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^ . Abstract 
» 

Analyzing the order of unobservable impulse in descriptor system leads to a new testing 
criterion for impulse observability, both the statement and the proof of which use only the 
J> ' original system data. 
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1 Introduction 

We consider the impulse observability of descriptor linear system [1], [2] 

Exit) = Ax{t) ^^^^ 
yit) = Cxit) 

where E, A & M"^", C G M™>^"-. Matrix E is singular, but matrix pencil sE — A is regular, i.e., 
(\et{sE—A) is a nonzero polynomial on s G C [3]. Comparing with standard linear system {E = I, the 
identity matrix), descriptor one is featured by having impulse behavior. The underlying mathematics 
for this phenomenon is that the following initial value problem of differential-algebraic equation 

Ex{t) = Ax{t), t>0; x{0)=w (1.2) 
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has no solution generally in the sense of classical differentiable function, and a generalized solution, as 
a mathematical model of the state response of the system to initial value, in the sense of distribution 
has to be adopted. It interprets the impulse behavior that the distributional solution may contain 
a linear combination of Dirac delta distribution S{t) and its distributional derivatives S^''\t), k — 

1,2, This linear combination is called impulsive term in the state response for convenience in 

following. For details, see [2, pp. 16-22] and the references therein. More recent works about 
distributional solution are [4] -[8]. 

System (1.1) is called impulse observable (see, e.g., [2, p. 43], [9]-[13], [6]), if for arbitrarily given 
but unknown initial value a:(0), the impulsive term in state response x{t), t > can be uniquely 
determined out from the measured, therefore known, output response y{t), t >0. 

A subtle point for impulse analysis in descriptor system is that explicit expressions of the impulsive 
terms in state and output responses directly using the system data, i.e., {E,A,C} for system (1.1) 
discussed here, do not exist. As a result, the existing impulse observability criteria, mainly collected 
in [2, Theorem 2-3.4, p. 43 ] and the dual form of [1, Theorem 4.2, p. 20], and/or their proofs 
have to employ indirect data which comes from some transformations to the original system, e.g., 
slow-fast decomposition based on Weierstrass canonical form, descriptor system structure algorithm 
[1], etc. Besides causing numerical difficulties, to deal with which some condensed forms based on 
numerically reliable orthogonal matrix transformations are developed [14]-[19], such type of answer 
to a control theory problem can not satisfy a theoretical interest as well. 

In this note, we are interested in obtaining an intrinsic impulse observability criterion in the 
sense that both its statement and its proof use the original system data directly. This concern is also 
motivated by seeking for a criterion with clear dynamical interpretation. The obtained new criterion 
is Theorem 2.4, the dynamical interpretation of which is the nonexistence of unobservable impulse 
of specific order, see Theorem 2.3. Main existing criteria for impulse observability in literature are 
special cases of our new one, see Corollaries 3.1 and 3.2. 



2 Main Results 



The frequency domain form of (1.1) with initial state x{0) = w can be written as 



sE-A 
C 



X{s)^ 



Ew 
Y(s) 



(2.1) 
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where X, Y denote the Laplace transforms of x, y respectively. 

We use notations M.'^(s) and R"[s] to denote the set of n-dimensional rational fraction vectors 
and the set of n-dimensional polynomial vectors respectively. Note that any rational fraction can 
be uniquely expressed as the sum of a strictly proper fraction and a polynomial. In following for 
F{s) e M.^{s) we always use Fa{s) and -Fp(s) to denote the strictly proper fraction pcirt 
and polynomial part respectively for notation convenience. 

For clarity we use the following Definition [7] . 

Definition 2.1 Let Vi e R", z = 0, 1, . . . ,p with Vp ^ 0. Then the distribution 

X{t) = 5{t)vo + 5^'\t)vi + ■■■ + 5^\t)vp 

is called an n-dimensional impulse of order p, denoted by deg(A) = p. The polynomial vector Yl^=o 
e in complex variable s of degree p, which is Laplace transform of impulse, is called frequency 

domain form of impulse, or simply, impulse. 

The Laplace inverse transform of a strictly proper fraction is a usual smooth function (a com- 
bination of exponential, triangular and polynomial functions of time variable t. see [20]). Therefore 
impulse observability of the system (1.1) means, in frequency domain language, that the polynomial 
part Xp[s) in X{s) can be uniquely determined out from Y{s). We write this fact as definition of 
impulse observability for clarity, although "common" definition is not so stated [2, p. 43], [?], [9], 
[10], [13], [6]. 

Definition 2.2 The system (1-1) is impulse observable, if the following equation 



sE-A 




Ew 




X{s) = 




C 








(2.2) 



on {w, X{s)) G M" X R"'(s) has no solution with Xp{s) nonzero. 

Remark 2.1 In other words, from Y(s) — CX(s) — one can conclude that X{s) does not contain 
polynomial part (frequency domain form of impulse, see Definition 2.1). 

Remark 2.2 We adopt the equation viewpoint. Here a dehberate point is that both vu and X{s) 
are seen undetermined simultaneously. It is this insight that results in an equal treating with p_i 
and po, - ■ ■ ,Pr in (2.12) and the introducing of the matrix (2.8). 
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Theorem 2.1 The system (1-1) is not impulse observable, if and only if the following equation 



sE-A 




Ev 




P{s) = 




C 








(2.3) 



on {v, P{s)) e M"- X M"'[s] has a solution with P{s) nonzero. 

Note the difference of and W[s]. 

Proof. The sufficiency is obvious by Definition 2.2 since P{s) G ]R"'[s] C ^""(5). 

Necessity. By Definition 2.2, there exists {w, X{s)) e M" x M'^(s) with Xp{s) nonzero such that 

sE-A 
C 

Firstly, = CXa(s) + CXp(s) e W^{s) i 

Secondly, the limit lim5_>oo sXx{s) exists 
Therefore 

Ew = {sE - A){Xa{s) + X^{s) 

= [{sE - A)Xa{s) - Eq] + [Eq + {sE - A)Xp{s)] 

forms a decomposition of strictly proper fraction plus polynomial. The uniqueness of such decompo- 
sition imphes Ew — [Eq + {sE — A)Xp(s)], which is equivalent to 

{sE-A)Xp{s) = E{w-q). (2.5) 

It follows from (2.4) and (2.5) that {w - q, Xp{s)) e x W[s] is a solution of (2.3). ■ 
We introduce the following technical notion. 

Definition 2.3 Let {v, P{s)) e R" x M.^[s] with P{s) nonzero and deg(P(s)) — r be a solution of 
the equation (2.3). Then v is called an unobservable impulsive initial state of the system (1.1), and 
P{s) is called an unobservable impulse of order r of the system (1.1). 

Lemma 2.1 The system (1.1) is impulse observable, if and only if it has no unobservable impulse 
of order < n — 1. 



{Xj,{s) + Xp{s) ^ 



mplies 



Ew 




CXp{s) = 0. 



(2.4) 



, denoted by q, and moreover \ims^(x>{sE — A)Xa{s) — Eq. 
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Proof. It follows from the regularity of the pencil sE — A that the solution of (2.3), if exists, will 
be of deg(P(s)) < n — 1. ■ 

Theorem 2.2 Let {v, P{s)) eW x W[s\ be a solution of (2.3) with deg(P(s)) ^r>l and write 

P{s) ^po-spi + --- + i-sYpr. (2.6) 

Then (pr-i, Pr) € M** x W^[s\ is a solution of (2.3) as well, where Pr G ]R"[s] is of deg{pr) — as a 
polynomial vector. 

Proof. Combining (2.3) and (2.6), we have 

r 

E 



i=0 



s'+^E - s'A 




Ev 




{-lypi = 




s'C 








Differentiating r times gives 

r\E-0 


and further 



-ly-'Pr-i + 



(r + ly.sE - r\A 
r\C 


















sE-A 




Epr-1 




Pr = 




C 








(2.7) 

I 

Note that (r + l)sEpr — imphes Epr = 0. ■ 

Prom an unobservable impulse P{s) of order r to initial value v, Theorem 2.2 constructs an 
unobservable one pr of order zero to initial value Pr-i- 

Form the following partitioned matrix 

E A 

E ■-. 



Ok{E,A,C)^ 



A 
E 



C 



(2.8) 



C 

for k — 2, . . . ,n + 1, where the other blocks not appearing are zero. It has k block columns and 
2k — 1 block rows, and then has the size {kn + {k — l)m) x kn. 



Theorem 2.3 The system (1-1) has no unobservable impulse of order < r if and only if 

ieink{Or+2{E, A, C)) = n{r + 1) + rank(£;). (2.9) 

Proof. Let deg(P(s)) < r and write P{s) — po — spi H h {—sYpr- Then (2.3) is equivalent to the 

following (r + 2) + (r + 1) = 2r + 3 equations 

Epi + Api+i = 0, i = -1, 0, . . . , r, (2.10) 

Cpi^O, i = 0,...,r (2.11) 

on {p-i,po, ■ ■ ■ , Pr) e M" X • • • X M", where we denote p_i = v and Pr+i = for notation convenience. 
The group of equations (2.10) and (2.11) can be rewritten into the following matrix form 



Or+2{E,A,C) 



P-1 

Po 

Pr 



(2.12) 



Therefore the system (1.1) has no r order unobservable impulse if and only if all solutions of (2.12) 
satisfy pi = 0, i = 0,1, ... ,r, i.e., 



■.Or+2{E,A,C) 



( 


P-1 




P-1 


> 






P-1 


< 


Po 


■.Or+2{E,A,C) 


Po 


= 


> = < 









Pr 




Pr 














> 


P-1 











= > 





> 



It is easy to see that 



Or+2{E,A,C) 



P-1 






= Ep-i. 



Then by computing dimensions of solution spaces of linear equations, (2.13) gives 

n(r + 2) - rank(C>^+2(£^, A,C))^n- rank(E) 
and the result follows immediately. ■ 
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Theorem 2.4 The following statements are equivalent: 

1) . The system (1-1) is impulse observable; 

2) . rank((9r+2(-E, A, C)) = n{r + 1) + rank(£') for some one r G {0, 1, ■ ■ ■ , n — 1}; 

3) . rsink{Or+2{E , A,C)) = n{r + 1) + rank(£;) for each one r e {0, 1, • • • , n - 1}. 

Proof. Follows from Lemma 2.1, Theorems 2.2 and 2.3. ■ 

Remcirk 2.3 Up till now, all statements and proofs use only the original system data, not involved 
in any transformation to the system. 



3 Implying Existing Results 

First, taking r — 0, the condition (2.9) gives the following consequence. 

Corollary 3.1 The system (1-1) is impulse observable, if and only if r&,vik{02{E, A,C)) 
rank(£'), i.e., 

E A 
rank E 
C 



n + rank(£'). 



n + 



(3.i; 



Remark 3.1 The condition (3.1) is a well known criterion for impulse observability [2, Eq. (2-3.6), 

p. 44], which is featured by using the original system data. To our knowledge, it is the only criterion 
of this feature in literature. Existing proofs (see [2, p. 44], [13], etc.) rely on some decompositions. 
Now we see that it guarantees the nonexistence of unobservable impulse of zero order. In quantitative 
aspect, equation 

E A 





E 
C 



P-i 



provides total information about unobservable impulsive initial states of zero order. 
Now we consider the Weierstrass canonical decomposition 



T 




sE-A 


I 




C 



s 



sL 



A, 









sN -L 



n2 



(3.2) 



7 



where T, S e R"^" are invertible, rii + n2 — n and A'" is nilpotent with index h (i.e., 7^ 0, but 

Nf" = 0). 

Lemma 3.1 The condition (2.9) holds if and only if 



rank(0^+2(A^, 4^, C2)) = n2(r + 1) + rank(Ar). 



Proof. Straightforward. 



Lemma 3.2 



rank(e)^+2(A^, C'2)) = n2(r + 1) + rank 



C2N 

hi 



CoN^ 



Proof. Through a series of elementary row and column transformations, matrix Or+2{N, ^n2jC2) 
can be transformed into 

In, 

■• 

jYr+2 

C2N 



^n2 





and the result follows. 



Lemma 3.3 The condition (2.9) holds if and only if 



rank 



C2N 

'2^ 



CoN^ 



rank(A^). 



(3.3) 



Proof. Follows from Lemmas 3.1 and 3.2. 
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Corollary 3.2 The following statements are equivalent: 

1) . The system (1-1) is impulse observable; 

2) . Condition (3.3) holds for some one r G {0, 1, ... ,n — 1}; 

3) . Condition (3.3) holds for each one r e {0, 1, . . . , n — 1}. 

Proof. Follows from Theorem 2.4 and Lemma 3.3 immediately. ■ 

Remcirk 3.2 When r > h — 2 with h the nilpotency index of N, the rank criterion (3.3) gives the 
condition [2, Theorem 2-3.4, (iv)]; When r = 0, it gives another condition [2, Theorem 2-3.4, (v)]. 
We see that these two well known criteria are only two boundary cases of Corollary 3.2. 
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